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Abstract. - The interplay of topological constraints, excluded volume interactions, persis- 
tence length and dynamical entanglement length in solutions and melts of linear chains and 
ring polymers is investigated by means of kinetic Monte Carlo simulations of a three dimen- 
sional lattice model. In unknotted and unconcatenated rings, topological constraints manifest 
themselves in the static properties above a typical length scale d t ~ (<f) being the volume 

fraction, I the mean bond length). 

Although one might expect that the same topological length will play a role in the dynamics 
of entangled polymers, we show that this is not the case. Instead, a different intrinsic length d e , 
which scales like excluded volume blob size f , governs the scaling of the dynamical properties of 
both linear chains and rings. In contrast to dt, d e has a strong dependence on the chain stiffness. 
The latter property enables us to study the full crossover scaling in dynamical properties, up 
to strongly entangled polymers. In agreement with experiment the scaling functions of both 
architectures are found to be very similar. 



Introduction. Topological constraints (TC) are commonly invoked in polymer theory to 
explain the dynamical behaviour of large molecular weight chains. In the reptation theory 
the long range TC are assumed to be responsible on a local scale for the conjectured 
anisotropic motion of chain segments, which are not allowed to cross neighboring chains. 
Obviously, in systems of linear chains this non-crossing constraint does not provide TC in 
a strict mathematical sense [^j. Strict topological constraints, however, do exist in systems 
of unconcatenated and unknotted ring polymers, where they strongly influence the static 
properties [^|-j8|. In this case, the TC are created by the non-crossing requirement and, as 
has been shown recently Q, a mass independent length scale dt- It is, hence, natural to 
ask whether dt provides the conjectured length scale, i.e. the "tube diameter" d e Q], for the 
dynamics of linear chains and rings. A link between the two concepts was indeed suggested by 
the fathers of the reptation concept, de Gennes (lj and Edwards 0]. This idea prompted the 
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Fig. 1 Fig. 2 

Fig. 1 - Ring size R = (Rq — Ro+N/2) 1 ^ 2 versus chain length N for dilute ((f) — > 0) and molten 
rings. The dilute rings are characterized by the same Flory exponent v = vq m 0.59 as their linear 
counterparts (not included). In the dense limit short rings form Gaussian chains of blobs [y = 1/2: 
solid line). With chain overlap, i.e. with increasing molecular mass, volume fraction and (most 
importantly) stiffness, the rings become more and more compact. 

Fig. 2 - Diffusion coefficient D(N, <f>, a) versus chain length N. We consider flexible (a — 0) and 
semiflexible (er = 3) linear chains (full symbols) and rings (empty symbols) in the dilute and the melt 
limit. In the dilute limit both architectures have the same Rouse-like dynamics (upper bold line). 
Contrasting to this the dynamics at high overlap depends strongly on both architecture and stiffness. 
Effective power law exponents D oc are indicated in the figure (dashed lines); (3 increases with 

N,(j> and a. Rings are always faster than corresponding linear chains and have smaller apparent 
exponents. The (3 « 2-slope — consistent with the reptation prediction jl|,|2) — is observed over more 
than one order of magnitude for semiflexible (a = 3) linear chains (lower bold line). 



present computational study, which extends our previous investigation on semiflexible rings 
on linear polymers. Concentrating on the scaling analysis of global quantities such as the 
diffusion constant and the radius of gyration, we demonstrate that the static and dynamic 
problems involve different length scales, with different dependence on chain stiffness and on 
density. The conjecture d e « d t is therefore incorrect, and the simpler proposal d e « £ @,@], 
where £ is the correlation blob size, accounts much better for the scaling behaviour. 

Computational method and raw data. As in our previous studies we use Monte-Carlo 
simulations of an extensively studied lattice model, the bond fluctuation model This 
model, in which the dynamics proceeds by local jump attempts, is thought to describe well 
the polymer dynamics on Langevin equation level |ll|] . In order to tune the persistence length, 
we impose a simple intramolecular potential which favors straight bond angles E(9) — a cos(#) 
where 9 denotes the complementary angle between two successive bonds. The stiffness a (in 
the range a = 0, 1, 2 and 3) is an energy scale, expressed here in units of the thermal energy 

In fig. I and fig. H we present two sets of raw data from our simulations — the ring size 
R and the self-diffusion constant D for linear chains and rings as functions of their molecular 
mass N. We have sampled a broad range of volume fractions (of occupied lattice sites) 4> from 
the dilute up to the melt [|ll| limit (</> ~ 1/2) with chain masses up to N = 1024 in order to 
vary the chain overlap sufficiently to test the proposed scaling scenarios. 



M. MULLER et al: TOPOLOGICAL CONSTRAINTS VS. ENTANGLEMENT LENGTH 



3 



In the bond fluctuation model the bond length I is not, strictly speaking, a constant, since 
a bond can be represented by lattice vectors of different lengths. The average bond length 
I, however, depends very weakly on stiffness and density, and will be considered as constant. 
The data for dilute chains can be used to define the effective or statistical segment length 
60 (er), obtained from fitting the radius of gyration of asymptotically long dilute linear chains 
with Rgyr = ba(a)N where vq ss 3/5 is the swollen chain exponent. &o(c) increases with 
stiffness up to 50%. Similarly, linear chains above the overlap concentration are Gaussian, 
and obey the scaling R gyr = &(</>, o^iV 1 / 2 . Local dynamics is characterized by a monomer 
mobility m((f>, a), obtained from the mean-square displacements (MSD) of center monomers 
at short times Again, this quantity turns out to be relatively insensitive to stiffness. In 
the dilute limit it may be equally be obtained from the center-of-mass motion. As shown in 
fig. H we find D = mo(a)/A for both architectures, with a mobility too roughly proportional 
to the acceptance rate. At higher densities we find m(4>,o~) « mo(a)m(4>) with to(0) = 1, i.e. 
density and rigidity effects decouple. 

While static and dynamic properties of both architectures resemble each other in the dilute 
limit and at moderate overlap, where both linear chains and rings form Gaussian chains of 
blobs - bold line in fig. [l| - marked differences appear at high chain overlap, as seen in figures 
[j] and ||. The local slopes in these figures, v = dlogR/dlogN and (3 = —dlogD/dlogN 
depend continuously on N, indicating that our data extensively cover the crossover regions. 
With increased chain overlap the rings become more compact (1/v —* 3) and their dynamics 
more entangled-like, i.e. the diffusion exponent f3 approaches 2 [Hfl . The stiffness dependence 
is very strong, especially for the diffusion constant. That the reptation prediction (3 = 2 for 
linear chains is indeed seen so clearly, is one of the important results of this work; this is 
in agreement with a very recent MD simulation p3| . Note that rings are always faster than 
their linear counterparts. 

Obviously, the slopes are only effective exponents which characterize the crossover between 
different regimes and do not capture asymptotic behaviour (e.g., the two slopes indicated in 
fig. [I] would then intercept which is unphysical) . The task is then to understand the effects on 
both architectures of increased chain overlap — specifically, the strong stiffness dependence 
— and to elucidate the underlying crossover scaling. 

Two intrinsic length scales and scaling attempts. There are two ]l5| ] natural intrinsic 
length scales — both chain length independent — which have been used to characterize the 
interactions in strongly overlapping polymer solutions and melts the 'diameter' p(<p) — 
l/VUt> and the 'blob' size £(<j>,b ) » b N^° where N t w (& i V) -1/(3I ' 0_1) [ s the number 
of monomers contained in the blob. Note that in contrast to p the blob size £ oc £>o _1 3 is 
strongly stiffness dependent. The prefactors of £ and iVg can be estimated independently from 
the static structure factor or, equivalently, from 6(</>, er) — C/N^ 1 ^ 2 |^,^|. 

We demonstrate now that p determines dt (Fig. H) and that the polymer dynamics is 
indeed characterized by one length scale d e ~ £ (Fig. ^) . Note that the figures include data 
from very different densities and rigidities and that no shift or fit parameter has been used. 
It is also important to remark that data points may be found in the different parts of the 
scaling functions, irrespective of their density or stiffness. This shows that the true intrinsic 
parameter is neither tf> nor a, but the proposed scaling variable. 

First, we have verified that the traditional scaling versus N/N/: § works for semi- 
flexible linear chains (not shown). Surprisingly, it does not work for rings This is an 
indication of the influence of TC on static properties. In fig. ||, we show that a proper scaling 
of the radius of gyration for rings is obtained by using the scaling variables R/d t and N/N t , 
where d t ~ 2h/\fUf>. The number N t of monomers between topological obstacles N t used in 
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Fig. 3 - Data collapse of the reduced ring size R/dt versus N/N t for rings at three representative 
densities (<f> = 1/2: full symbols; <f) — 1/8: large, thin symbols; (f> — 1/32: small symbols) at different 
a as indicated in the figure. The successful scaling attempt supposes d t ~ 25p and self-consistently 
(full line on the left) Gaussian chains of blobs at small distances, i.e. N t ~ 4O5(p/&(0, cr)) 2 . The slope 
v = 1/4 for the predicted LA asymptotic behaviour for larger chains is indicated on the right (dashed) 
along with the apparent v — 1/3-slope (dashed-dotted) from fig. ji] which fits a broader range of data. 
Inset: The bond-bond correlation function (eo ■ e*o+ n ) shows striking anti-correlations for rings (bold 
line) at large overlap (iV = 1024, <f> = Q.5,a = 3) in contrast to linear chains (dashed). 



fig. H is obtained from dt ~ b((f>, a)^ 1 ^ 2 where we assume Gaussian statistics at distances 
smaller dt- This is motivated by fig. [j] and is consistent with the scaling (left side of fig. ||). 
From this successful scaling collapse we may conclude that the length generated by the TC 
is d t « 25/ This yields d t « 22 for = 0.5 and N t (a = 0) w 173 and N t (a = 3) « 63. 
This length is both mass and stiffness independent. The prefactor is somewhat arbitrary and 
depends on the asymptotic slope on the right which is numerically unclear (see below). 

We turn now to the scaling of the diffusion data, and demonstrate, that polymer dynamics 
is indeed characterized by one length scale d e which scales like the blob size £. In fig. [|, 
the mass has been rescaled by a number N e assumed to be proportional to the number of 
monomers per blob N%. The scale for the vertical axis is given by D e — m((f>,a)/N e , i.e. we 
assume Rouse-like behaviour for weakly overlapping chains. For small N all data collapse 
irrespective of architecture and stiffness on the asymptotic slope with (3=1. This shows that 
the definition of D e is self-consistent and also that the mobility correction m(4>, a) was chosen 
appropriately. For N 3> N e we find that the linear chain data follow the reptation prediction 
D/D e = (N e /N) 2 for nearly two orders in magnitude! From the known blob size and the 
intercept of both asymptotes one estimates N e /N^ w 15 for the linear chains. Note that our 
scaling is consistent with the traditional density crossover scaling assumption B and recent 



experiments on linear chains 114,151. Also, the present scaling, in contrast to that in Ref. |R 
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Fig. 4 - Scaling attempt for diffusion constants of linear chains and rings (shifted upwards for 
clarity) supposing d e ~ £. D/D e is plotted versus N/N e where we use D e — m(<f>,a)/N e and 
iV e = 32.8(&o 3 0) _1//< ' 3l ' o_1 - 1 . Data for four rigidities (as indicated in the figure) and three densities 
(as in fig. pi same symbols) have been included. The data collapse was achieved without any free 
parameter! The linear chain data follow the reptation prediction for nearly two orders of magnitude. 
The prefactor of N e was defined such that both linear chain asymptotes (bold and dashed lines) 
intercept at D/D e = N/N e = 1. 



does not involve any arbitrary shift parameters. 

The universal functions for the dynamical properties of both architectures are similar. In 
agreement with experiment jj, 16 ring dynamics becomes entangled- like at slightly higher 
chain length. This is no surprise and may be attributed to their smaller size P]. However, 
the asymptotic behaviour is less obvious. The best fit for the available data (3 1.7, may 
be regarded as lower bound. Motivated by recent experiments |l(| which show the same 
asymptotic behaviour for both architectures we have also fitted the ring data with (3 = 2 
yielding a reasonable agreement (upper dashed line in fig. ^). From the intercept of the 
unentangled and entangled asymptote one concludes that the ring crossover occurs roughly 
at 3 — 4 times larger chain length than in their linear counterparts — again in agreement 
with experiment |]l6f . This finally explains why the (3 exponents of rings at given <p and a are 
smaller. 

Interpretation of scalings results and new queries. Both presented scaling plots are consis- 
tent with two eminent mean-field descriptions in polymer physics — the lattice animal (LA) 
proposal for dense rings and the reptation concept for polymer dynamics. They pro- 

vide additional information in that they identify the phenomenological length scales implicit 
in both theories. 

Our compact rings at high chain overlap behave like LA in a self-consistent network im- 
posed by the TC of neighboring rings, i.e. they may be described on the same footing as 
isolated rings in a network of fixed topological obstacles of density l/d t 3 H. If the ring size 
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R becomes larger than d t the rings are forced to retrace their paths and the fractal dimension 
\jv becomes that of a strongly branched object. Within the single chain model it is known 
that the rings are characterized by the exponent v = 1/4 within an intermediate, but broad 
chain length window ||. It is this suggestion we have followed in fig. || to fix the prefactor of 
dt- Note that the rings presented in fig. |^ are strongly overlapping, with R much larger than 
the distance between chains (see (§), fig. 7), and that we refer to 'compactness' in the scaling 
sense, v < 1/3. 

The effective topological interactions can be visualized directly with the bond-bond vec- 
tor correlation function shown in the inset of fig.[|. This function has a pronounced anti- 
correlation. The negative minimum indicates that the ring has to fold back after n ~ 10 
monomers to retrace its path. The position and the depth of the dip increase with <f> and a, 
but are chain length independent which underlines the AMndependence of dt- The dt ~ p 
relation reflects the topological origin of the interaction. Indeed p characterizes the volume 
per chain segment which is invariant with regard to chain conformation, chain length and 
stiffness. The advantage in using stiff chains is related to the fact that dt does not depend on 
stiffness: semiflexible rings 'waste' fewer monomers on short length scales, may explore larger 
distances and become more compact. Tuning the stiffness provides therefore a much more 
efficient route to increase the chain overlap, required for the mean-field picture to hold, than 
a chain length variation [|| . 

Confirming older j^] and more recent simulations |^,|l3| we find that increasing the stiffness 
is an even more efficient trick for exploring the strongly entangled dynamics. This is now 
rationalized by the scaling d e ~ £ in fig. ^ which implies that the strong stiffness dependence 
of the number of entanglements per chain N/N e oc &o 3 ' 9 0]- We emphasize, however, that we 
have not demonstrated in this study the reptation proposal. The j3 — 2 asymptote of fig. ^ is 
not sufficient to claim anisotropic motion of the chains along their contour || which requires, 
e.g., a more detailed analysis of the MSD. This must take into account the non-asymptotic 
behaviour in the MSD generated by tube renewal [pl| and constraint release effects. The 
latter mechanism is likely to be the dominant relaxation mode for our entangled rings. What 
has been demonstrated here in a computational study covering, for the first time, several 
orders of magnitude is (i) the existence of two distinct dynamical regimes, both in linear 
chains (consistent with j|, |^, O, ) and ring polymers, (ii) that the dynamical crossover is 
characterized by one length scale only, (iii) that this length scale is the size of the excluded 
volume blob, i.e. a properly defined static quantity |l9| , (iv) that, surprisingly, the crossover 
scaling does not depend on the architecture. However, the finding d t 3> d e for the strongly 
overlapping chains in our simulations may explain naturally why the length scale dt appears 
to play no role in the dynamics of rings of experimental |lq| and computational relevance. 

Summary. In this letter, we have elucidated, by means of a numerical experiment, the 
scaling of static and dynamic properties of strongly entangled linear chains and rings with 
regard to density and stiffness variation. Our results indicate that topological constraints 
manifest themselves differently in the statics and in the dynamics. The topological length 
which governs the crossover in the statics of rings is essentially independent of stiffness. This 
contrasts with the strong stiffness dependence of the dynamical entanglement length, which 
governs the crossover to rcptation-likc behaviour. Although the high molecular weight dif- 
fusion exponent agrees with the reptation prediction, we have not presented rigorous proofs 
for the anisotropic motion of polymer chains. A very puzzling question in this context is the 
similarity between the scaling functions for linear chain and ring dynamics which is the focus 
of forthcoming work. 
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